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In this paper we study linear and non-linear Fourier-Pade approximation 
for Angelesco systems of functions. This construction is similar to that of 
Hermite-Pade approximation. Instead of considering power series expan- 
sions of the functions in the system, we take their expansion in a series of 
orthogonal polynomials. 

In [£] and 0, S. P. Suetin obtained convergence results for rows of Fourier- 
Pade approximation extending to this setting classical results of the theory 
of Pade approximation. 

Diagonal sequences of Fourier-Pade approximation were studied by A. 
A. Gonchar, E. A. Rakhmanov, and S. P. Suetin in ,2J when the function 
to be approximated is of Markov type; that is, the Cauchy transform of 
a measure supported on the real line. They give the rate of convergence 
of diagonal sequences of linear and non-linear Pade approximants in terms 
of the equilibrium measures of a related potential theoretic problem. We 
generalize those results to the case when a system of Markov functions is 
given defined by measures whose supports do not intersect. 

Let A4(A) denote the class of all finite, Borel measures with compact 
support consisting of an infinite set of points contained in an interval A of 
the real line K. Given a G M(A), let 



be the associated Markov function. Let A&, k = 1, . . . , m, be intervals of the 
real line such that 



and o~k S At(Afc), k = 1, . . . , m. We say that a = (a%, . . . , a m ) forms an 
Angelesco system of measures and (a±, . . . ,a m ) is the associated Angelesco 
system of functions. 
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Let erg € A^(Ao). Likewise, we will assume that 

A o nA fc = 0, k = l,...,m. 

Consider the sequence {£ n }, n G Z + = {0, 1, 2, . . .}, of orthonormal polyno- 
mials with respect to <jq with positive leading coefficient. Take a multi-index 
n = (m, . . . , n m ) € Z™. Set 

|n| = m -\ h n m . 

Let Q n , P n ,i, • • • , Pn,m, be polynomials such that: 

i) degQ n < \n\,Q n 0,degP nJ < |n| - l,j = l,...,m. 

ii) For each j = 1, . . . , m , and k = 0, . . . , |n| + rij — 1 



The |n| + 1 coefficients of Q n satisfy a homogeneous linear system of |nj 
equations given by 

CkiQnVj) = 0, j = l,...,m, k = |n|, . . . , |n| + rij - 1 . 

Therefore, a non-trivial solution is guaranteed. 

In Section |^1 we will prove that every solution to i)-ii) has deg Q n = 
|n|. This being the case, (Q n) Pn,i, ■ ■ ■ , Pa,m) is uniquely determined up 
to a constant factor. In fact, let us assume that (Q n ,Pn,i, ■ ■ ■ ,Pn,m), and 
(Qn, Pn,i, ■ ■ ■ , Pn,m), are solutions of i)-ii). Without loss of generality, we can 
assume that Q n and Q n are monic (with leading coefficient equal to one). 
Obviously, if Q n - Q n ^ then (Q n - Q n , P n 

,1 Pn,li • • • ) Pn,m Pn,m) IS 

also a solution with degQn — Qn < |n| which contradicts our assumption. 
Hence Q n = Q n and by Q it follows that P n j = P n j,j = 1 . . . , m . 



The rational vector function f . . . , -^-J constructed from any so- 
lution of i)-ii) is called the n-th linear Fourier-Pade approximant for the 
Angelesco system (di, . . . ,a m ) with respect to uq. We shall see that for 
all n £ Zip the linear Fourier-Pade approximant of an Angelesco system is 
unique. 

Non- linear Fourier-Pade approximants are determined as follows. Given 
n € Z™, we must find polynomials T n , 5 nj i, . . . , S njm such that 

i') degT n < |n|,T n ^ 0,deg(5 nj ) < |n| - 1, j = 1, . . . ,m . 
ii') For each j = 1, . . . , m , and k = 0, . . . , |n| + rij — 1 




Notice that 



n|-l 



(1) 
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This system of equations is non-linear in the coefficients of the polynomials. 
We shall prove that for each n E Z!p, the system has a solution but we have 
not been able to show that it is unique. For any solution of i')-ii'), the vector 

rational function (^-^-, ■ ■ ■ , S ^' m ^j is called an n-th non- linear Fourier-Pade 

approximant for the Angelesco system (<7i, . . . , a m ) with respect to co- 
in this paper, we obtain the rate of convergence (divergence) of linear and 
non-linear Fourier-Pade approximants for Angelesco systems such that the 
measures ctq, . . . , a m are in the class Reg of regular measures. For different 
equivalent forms of defining regular measures see sections 3.1 to 3.3 in jSJ. 
In particular, <7o € Reg if and only if 

lim \i n (z)\ l/n = exp{g n (z; oo)} , 

n 

uniformly on compact subsets of the complement of the smallest interval 
containing the support, supp(cro), of uq and gn (-) oo) denotes Green's func- 
tion for the region = C \ supp(<7o) with singularity at oo. Analogously, 
one defines regularity for the other measures a±, . . . , a m . In the sequel, we 
write (do; o"i, . . . , a m ) G Reg to mean that £ Reg, k = 0, . . . , m . The 
system (a±, . . . , a m ) will be used to construct the Angelesco system of func- 
tions whereas do will determine the system of orthogonal polynomials with 
respect to which the Fourier expansions will be taken. Therefore, for all 
< j,k < m , we assume that 

Aj n A fc = , j^k. 

In Theorems 1 and 2 below, we find the rate of convergence of the |n|th 
root of the error of approximation of the functions Uk by linear and non- 
linear Fourier-Pade approximants, respectively. The answers are given in 
terms of extremal solutions of certain vector valued equilibrium problems 
for the logarithmic potential. Before stating Theorems 1 and 2, we need to 
introduce some notation and results from potential theory. 

Let Ffc, k = 1, . . . , N, be (not necessarily distinct) closed bounded intervals 
of the real line and C = (cj fe) be a real, positive definite, symmetric matrix of 
order N. C will be called the interaction matrix. By M\{F^), k = 1, . . . , N, 
we denote the subclass of probability measures of M.{F}.) and 

Mi = Mi(Fx) x • • • x Mi(F N ) . 

Given a vector measure /i £ A4i and j = 1, . . . , N, we define the combined 
potential 

N 

(2) Wf{x) = Cj,kV^ (x) , x € Aj , 

k=l 

where 

W*(x) = f log Tadjik (*) 
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denotes the standard logarithmic potential of fik- We denote 

w% = inf{W^(x) : x G Aj} . 

In Chapter 5 of [3| (see Propositions 4.5, 4.6, and Theorem 4.1) the au- 
thors prove (we state the result in a form convenient for our purpose) 

Lemma 1. Let C be a real, positive definite, symmetric matrix of order N. 
If there exists ~p, = (p l5 . . . ,~Pn) G M.\ such that for each j = 1, . . . , N 



x G supp(/j,A , 



then \i is unique. Moreover, if Cj t j~ > when Fj n Fj ^ then fx exists. 

The vector measure Jl G M\ is called the equilibrium solution for the 
vector potential problem determined by C on the system of intervals Fj , j = 
1,...,N. 

In the sequel A = A(pi, . . . ,p m ) C 
multi-indices such that 



lim —^-r 
neA n 



is an infinite system of distinct 
j = 1, . . . ,m. 



Let us define the block matrix 

Ci 

where 



Cl,2 
^2,2 



( 2p? 
P2P1 



P1P2 
2p| 



P2Pm 



2?& 



/ 



-p m (l +Pm)} , 



\ PmPl PmP2 ■ 

and £1,2, £2,1, £2,2 are diagonal matrices given by 

Ci,2 = C 2 ,i = diag{-pi(l +pi), -P2(l +P2), ' 
and 

C 2 , 2 = diag{2(l + Pl ) 2 , 2(1 + p 2 )\ ■ ■ ■ , 2(1 + p m ) 2 } . 
C\ satisfies all the assumptions of Lemma ^ on the system of intervals Fj 



A 



i,3 



1. 



, m, F i 



A 



0,3 



m + 1, . . . , 2m, including Cjk > when 



Fj nFj ^ 0. The only non-trivial property is its positive definiteness and we 
shall prove this in Section 2. Let ~p = ~p(Ci) be the equilibrium solution for 
the corresponding vector potential problem. We have 



Theorem 1 . Let (00 J 0°i , ■ ■ 

indices A = A(pi, . . . ,p m ). 



Let 



G Reg and consider the sequence of multi- 
Pn ' m ^ , n G A, be the associated 



sequence of linear Fourier-Pade approximants for the Angelesco system of 
functions (d\, . . . , d m ) with respect to gq. Then, 

= Gj{z), j = l,...,m, 



(3) 



lim 

neA 



P nj(z) 



Qn(*) 
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uniformly on each compact subset ofC \ (U^L Afc), where 

f(z)-uj)/ Pj 
and the combined potentials are defined by 



exp[(Whz) 



Set 



using C\ 



Gf = {xeC\(uT =0 A k ): ±[ur : 



W?(x)) > 0} 



An immediate consequence of Theorem ^ is 
Corollary 1. Under the assumptions of Theorem^ 

dj , 7 = 1 



p . 

lim 



m . 



ISA Qr 

uniformly on compact subsets of Gj and diverges to infinity at each point 
ofGj. 

Non-linear Fourier Pade approximants require the solution of a different 
vector potential equilibrium problem. Let 



Co 



c 2 



2.2 



where Ci,i, (+^2, £2,1 are as before and 

2m(l+pi) 2 



r 2 

u 2,2 



V 



m+1 
-2(l+p 2 )(l+pi) 
m+1 



-2(l+p m )(l+pi) 
m+1 



2(l+ Pl )(l+p 2 ) 
m+1 
2m(l+p 2 ) 2 
m+1 



-2(l+p m )(l+p 2 ) 
m+1 



-2(l+pi)(l+p m ) \ 

m+1 
-2(l+p 2 )(l+p m ) 

m+1 



2m(l+p m ) 2 
m+1 



/ 



C2 is a real, positive definite, symmetric matrix of order 2m. We take the 
system of intervals Fj = Aj,j = 1, . . . , m, Fj = Aq, j = m + 1, . . . , 2m. C2 
does not satisfy that Cj ^ > when Fj n Fj 7^ 0. In Theorem ^ of Section 
3, we prove that the corresponding equilibrium problem has at least one 
solution and that C2 is positive definite. Therefore, according to Lemma ^ 
the solution is unique. Let Ji = JLiC?) be the equilibrium solution for the 
corresponding vector potential problem. In Lemma El we show that for each 
n € there exists at least one non-linear Fourier- Pade approximant but 
we have not been able to prove that it is unique. We have 



Theorem 2. Let (<to; o~x, • • • , o~ m ) G Reg and consider the sequence of multi- 

Sn^l Sri,? 



indices A = A(p 



,Pr, 



Let 



, n G A, be an associated 



sequence of non-linear Fourier-Pade approximants for the Angelesco system 
of functions (ax, ■ ■ ■ ,a m ) with respect to do. Then, 

l/H 

= Hj(z), j = l,...,m, 



(4) 



lim 

neA 



adz) 



S n,j(z) 



Tn(z) 
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uniformly on each compact subset ofC \ (U™ =0 Aj), where 
H j (z) = exp((Wf(z)-cJf)/p j ) , 

fx = ~fi(C2 ) and the combined potentials Wj are defined by @ using C2 ■ 

Notice that the limit only depends on A and not on the non-linear Fourier- 
Pade approximants selected (in case that they were not uniquely deter- 
mined). Set 

Hf = {x G C \ (U™ oA, : ± (Jf - Wf(x)) > 0}. 
As a consequence of Theorem |2 we obtain 
Corollary 2. Under the assumptions of Theorem^ 

lim^ = ff i , j = l,...,m, 

neA J n 

uniformly on compact subsets of and diverges to infinity at each point 
ofHj. 

Section |2 is dedicated to the proof of Theorem^ and Section |3] to that of 
Theorem^ Section0]is dedicated to the justification of LemmaE]as stated 
here since in |3] the assumption Cj k > if Fj n Fj. ^ is assumed in general. 
In the sequel, we maintain the notation introduced above. 

2. Proof of Theorem [U 

From the definition of the linear Fourier-Pade approximant immediately 
follows that for each j = l,...,m 

(5) J x k (Q n (x)dj(z) - P n>j (x))da Q (x) = 0, k = 0, . . . , |n| + rij - 1 . 

Since the function Q n {z)dj(z) — P n j(z) is continuous on A , from © we 
have that Q n (z)dj(z) — P n j(z) has at least |n| +rij sign changes on Ao- 

Let W n j be the monic polynomial whose zeros are the points where 
Qn(z)&j(z) — Pnj(z) changes sign on the interval Ao- Obviously, deg W a j > 
|n| + nj and 

Qn(z)dj{z) - P n j{z) - . 

J G H(C \ supp(<Tj)), j = l m, 

Wn,j{Z) 

is analytic on the indicated region. Thus, linear Fourier-Pade approximants 
satisfy interpolation conditions on Ao- A similar statement holds for the 
non-linear Fourier-Pade approximants. In our proofs, we will use certain 
orthogonality relations satisfied by vector rational interpolants. 

Lemma 2. Let (ax,...,a m ) be an Angelesco system, n = (m,...,n m ) G 
Z™, and (u> n ,i> ■ • • , ^n,m) a system of polynomials such that degw n j > |n| + 
rij, j = 1, . . . , m whose zeros lie on an interval Aq, AqHAj = 0, j = 1, . . . , m. 
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Let (^-, • • • , ^f 2 ") ^ e a vector rational function such that degp nj j < |n| 
l,j = 1, . . . ,m,degg n < \jx\,q u ^ 0, and 

(6) J - -— £n(C\supp{aj)), j = l,...,m. 



Then 



(7) x kJHty^da j {x) = 0, k = 0,l,...,n j -l, j = l,. 

Consequently, degq n = |n| wii/i exactly rij simple zeros in the interior of 
Aj (in connection with intervals of the real line, the interior refers to the 
Euclidean topology of the real line) and degw n j = |n| + rij, j = 1, . . . ,m. 
Let Qn — QnjQnji where q n j is the monic polynomial whose zeros are those 
of q n lying in the interior of Aj . Then 

(8) 9,w - ^# = , ^. (z) , . f tMMsUfr) . 

1n{z) qij{z)q n ,j(z) J z-x w nJ (x) 

Proof. Notice that © and the assumption on the degrees of the polyno- 
mials q n ,Pn,j, and w n ,j imply that for j = 1, . . . ,m, and k = 0, . . . , rij — 1, 

W n ,j{z) \Z l ) 

Let Tj be a closed, smooth, Jordan curve that surrounds Aj such that all the 
intervals Aj, i ^ j, i = 0, . . . ,m, lie in the unbounded connected component 
of the complement of Tj. By Cauchy's Theorem, Cauchy's Integral Formula 
and Fubini's Theorem, it follows that 

Q= J_ f z kQn(z)d j (z)-Pnj{z) dz = 

2m J r . Wn,j{z) 
1 f zk qn (z)aj(z) ^ _ J_ f zk PnM dz = 

' J n,j(Z) 



2m J r Wn,j(z) Jr. w 



— I / ^ M dz= fx k ^Ldo 3 {x), 

2m J T . w n j(z) J z-x J w n:j (x) 

for k = 0, 1, . . . , nj — 1 and j = 1, . . . , m. Therefore, (Q) follows. 

Using standard arguments of orthogonality, from (JJJ) we obtain that q n 
must have at least rij sign changes in the interior of Aj and, consequently, 
at least rij zeros of odd multiplicity. Since degg n < |n|, we have that 
deg^ri = l n |) that all its zeros are simple and they are distributed in such a 
way that exactly rij lie in the interior of Aj. 

Assume that degw n j > |n| + nj for some j. Then 

w n,j(z) \z 2 
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This implies that ||2J) holds for all k = 0, . . . , rij. In turn, this means that q n 
has at least rij + 1 zeros in the interior of Aj against what was just proved. 
Therefore, degw n j = |n| + nj. 

Set q n = q n ,jQn,j, where q n j is the monic polynomial whose zeros are 
those of q n lying in the interior of A.,-. Notice that q n jdaj/w n> j is a real 
measure with constant sign on A.,-. For future reference, notice that with 
this notation the orthogonality relations may be expressed as 

(9) / x k q n:J (x)\q nJ (x)\ daj ^\ =0, k = 0, 1, . . . , rij - 1 . 

Hence, for each j = 1, ... ,m, q n j is the monic orthogonal polynomial of 
degree rij with respect to the varying measure |5°' J \ do~j . 



Notice that 



W n ,j(z) \z, 

Choose Tj as before. Using Cauchy's integral formula, Cauchy's Theorem, 
and Fubini's Theorem, we obtain that for each j = 1, . . . , m 

[qn,j(qn&j - Pn,j)](z) __ l f [q n ,j{q n&j Pn. 



w n j{z) 2m J r . w n ,j(()(z-() 

1 f ( «" )(C) -d C da J (x)= [ q ^M^ daj{x) 



2m J T . w n j(()(z -()((- x) J Z-XW n j(x) 

which is equivalent to l|8|), □ 
The vector rational function (^^, • • • , ^p 1 ) is called a multipoint vector 
Pade approximant of the Angelesco system (cti, . . . , a m ). According to Lem- 
ma|2]a necessary condition for their existence is that degu; n j < |n| +nj,j = 
1, ... ,m. Solving a homogeneous linear system of equations one sees that 
this condition is also sufficient. When degw n j = |n| + rij,j = l,...,m 
uniqueness follows because then degg n = |n| as we have seen. 

Remark 1. Applying this Lemma to linear Fourier-Pade approximants, we 
have that deg(Q n ) = |n|. Thus, for each n € Z+, they are uniquely deter- 
mined as claimed. 

Let us return to linear Fourier-Pade approximants. In this case, Q n = 

Qn, Qn,j = Qn,j,Qn,j = Qn,j and W n j = W n j . 



Lemma 3. For each j = 1, . . . , m, and k = 0, . . . , |n| + n 

V ' J \Qn,j(t)\ \J \t~x\ \W n ,j(x)\ 3K 

Moreover, degW n j = |n| + rij,j = 1, . . . ,m; that is, Q a (z)aj(z) — P n j(z) 
has exactly |n| + rij sign changes in the interior o/Aq. 
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Proof. From JHJ) and the definition of the linear Fourier-Pade approxi- 
mant, IjlflJI follows directly. The assertion concerning the degree of W n ,j is 
also contained in Lemma |2j □ 

Let {hi} C M(JC) be a sequence of measures, where K, is a compact subset 
of the complex plane and \i G M.(K). We write 

* lim hi = fi , /i£M (JC) , 
if for every continuous function / G C(JC) 



lim y /djtij = J fdfi ; 

that is, when the sequence of measures converges to fx in the weak star 
topology. Given a polynomial qi of degree I > 1, we denote the associated 
normalized zero counting measure by 

u n = 7 Z 5 * ' 

where (5^ is the Dirac measure with mass 1 at x (in the sum the zeros are 
repeated according to their multiplicity). 

In order to prove our main results we need Theorem 3.3.3 of 5 a . We 
present it in the form stated in which is more adequate for our purpose. 
In PP, it was proved under stronger assumptions on the measure. 

Lemma 4. Let {<f>i},l € T C Z+, be a sequence of positive continuous 
functions on a bounded closed interval A C M, a G Reg n A4(A), and let 
{qi},l G r, be a sequence of monic polynomials such that degqi = I and 

q l (t)t k cj )l (t)do-(t)=0, k = 0,...,l-l. 

Assume that 

lim — - log - — — - = v(x), 
lev 21 6 |0,(z)| v h 



uniformly on A. Then 



* lim v„, 



and 

lim f / |%| 2 z d/u J = e~ w , 
where V G A^i(Ai) is the equilibrium measure for the extremal problem 

V v {x) + v(x) 
in the presence of the external field v. 



= u, x G supp{v) , 
> lu, x G Ai , 



Using this result, we can obtain the asymptotic limit distribution of the 
zeros of the polynomials Q n j and W n j. 
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Theorem 3. Let a±, . . . , a m ) G Reg and consider the sequence of multi- 
indices A = A(j?i, . . . ,p m ). Then, for each j = 1, . . . , m 

* n m fQ»j = ~Pj > * i™ ^ n ,, = ~P m +j , 

nGA J neA J 

where Ji = Jl{C\) £ Mi is the vector equilibrium measures determined by the 
matrix C\ on the system of intervals Fj = Aj,j = 1, ... , m, Fj = Ao, j = 
m + 1, . . . , 2m . 

Proof. The unit ball in the cone of positive Borel measures is weakly 
compact; therefore, it is sufficient to show that the sequences of measures 
{ u Qn,j} anc i { u W n ,j}-, nGA, have only one accumulation point which coincide 
with the components of the vector measure ~fi(Ci) respectively. Let A' C A 
be a subsequence of multi- indices such that for each j = 1, . . . , m 

* lim vn . = Vn , * lim uw . = v m +i . 

neA' gnj 3 ' neA' Wn ' 3 +] 

(Notice that Uj G Ad\(Aj),j = l,...,m, and Vj G .Mi(Ao),j = m + 
1, . . . , 2m.) Therefore, 

(11) lim JQ„j(z)|^ = exp(-^(z)), 

nGA' 

uniformly on compact subsets of C \ Aj , and 



(12) lim \W n j(z)\^+^ = exp(-V Um +j(z)), 

neA' 

uniformly on compact subsets of C \ Ao- 

For each fixed j = 1 . . . , m, the polynomials Q n ,j satisfy the orthogonality 
relations ©. Using JTTJ) and (JT2J) it follows that 

lim — log = _ 1 + W ^ (x) + V (x), 

uniformly on Aj. By Lemma 0] za- is the unique equilibrium measure for the 
extremal problem 

(13) V ^(x) + J2^V^(x)-^^V^(x)>e j , x E Aj , 
with equality for all x G supp(z^). Additionally, 




i 

2n j 



(14) lim, [ / IQn.(x)| 2 ^^Pj = e->. 

On the other hand, for each fixed j = 1, ... ,m, the polynomials W n j 
satisfy the orthogonality relations (|1U|) and we can apply once more Lemma 
H Notice that for all t G A 

, s f \Ql,j( X )\ \Q n ,j(x)\d(Jj(x) 1 /' ^| 2 IQn,j(g)ldg-j(a) 

( J y i*-x| i^ nj (x)i |w nij (x)i ' 
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where Sj = inf{|i — x\ : t G Ao, x G Aj} and 
y2 



(16) 



1^0*01 |Qnj(x)|dt7 3 -(g) _1_ 

\t-x\ \W nJ (x)\ -6* 



\Q 



"J 1 



,2 \Qn,j( x )\ da j( x ) 
\Wnj(x)\ 



with <5* = max{|i - x\ : t € A , x G Aj}. From (|TT|). (JT5J), dHJ), (fT5J). and 
(|lfj|l. we obtain 

1 



lim - .. 
neA' 2(|n| + rij 



log 



lgnj(g)| 



.I' 



?n,j(t)| 2 |Q n ,j(t)Mo-j(t) 
VK n j(t) 



Pi 



Pi 



2(1 + Pj ) ' 1 • 

uniformly on Ao- Using Lemma 01 ^m+j is the unique extremal solution for 
the equilibrium problem 



(17) V Um +i{x) 



P.i 



-V"i(x) + 



P.I 



1+Pj 3 



> 



'm+j j 



x G A c 



2(1+ P,) 

with equality for all x G supp(v m +j). 

Rewriting (|13|) and (|17jl conveniently, we see that the vector measure 
(z/i, . . . , z^m) G M.\ is the unique solution for the vector equilibrium problem 
determined by the system of extremal problems 

(18) 2p)V v i (x) + ^PjPkV* (x) - Pi (l + Pj )V Vm +i (x) > w is x G Aj , 

(2pj9j = Uj) with equality for all x G supp(^), and 

(19) 2(1 + Pj ) 2 V v ™+i (x) - Pj (l + Pj )V v * (x) > co m+j , x G A , 

with equality for all x G supp(f m +j). That is, it is the equilibrium measure 
JZ G Mi for the vector potential problem determined by C\ on the system of 
intervals Fj = Aj, j = 1, . . . , m, Fj = Aq, j = m + 1, . . . , 2m. The condition 
Cj,k > if Fj n .Ffc ^ is fulfilled. According to Lemma ^ this equilibrium 
vector measure is uniquely determined if C\ is positive definite. Let us prove 
this. 

For j G {1, . . . , m} the principle minor of order j of Ci is 



detCC? 5 ) = (p 1 ---p i ) 



2 1 
1 2 

1 1 



(Pi-"Pj) 2 (j + l) >o. 



For j G {m + 1, . . . , 2m} the principle minor of order j of Ci can be 

calculated as follows. For each k = 1, . . . , m, factor out pk from the feth row 

(i) 

and A;th column of C\ . From the row and column m + k,k = 1, . . . ,j — m, 
factor out 1 + pt- In the resulting determinant, for each k = 1, . . . , j — m, 
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add the kth row to the (m + k)th row and then to the resulting determinant 
add the kth column to the (m + k)th column. We obtain 



det(c[ j) ) = [p! • ■ -p m (l + P1 ) • • • (1 + Pj - m )f 



2 1 
1 2 

1 1 



[Pi • • -Pm(l +Pl) ■ 
With this we conclude the proof. 
Proof of Theorem^ From l| 



[l+Pj- m )] 2 {m+l) >0. 



□ 



, l|15j). and (|16|). the asymptotic behavior 
of the function — ^ depends on the behavior of W^j, Qnj, Qn,j, arid 
7 n j, where 



mm ■ 



7i 



\Q(x)\ 



2 \Qn,j{x)\d(Tj(x) 
\Wnj{x)\ 



: Q(x)= x n * + 



\Qn,j(X i 



1 2 \Qn,j(x)\daj(x) 



From Theorem El for each j = 1, . . . , m, we have 



(20) lim |W ni ,-(x)| 1/|n| = exp{-(l +p.-)^(i)}, 

nGA ' 

uniformly on compact subsets of C \ Ao, and 
(21) 



lim [Q n ,i(x)r /|n| = exp{- Pj V^(x)}, 

neA 

uniformly on compact subsets of C \ A,-, where /Z = ~p(Ci) Using Q14|) (see 
also parenthesis after <|18|)V it follows that 

i/M 



(22) 



1 



lim 



exp{— 2pj9j} = exp{— Uj/pj} . 



Combining ©, l(2D|>. (GU, and (J22J), we conclude that © holds true 
uniformly on compact subsets of the indicated region. □ 

3. Proof of Theorem [2] 

We begin by proving the existence of non-linear Fourier-Pade approxi- 
mants. 

Lemma 5. Given (<7q; <ti, . . . , cr m ), for each n £ there exists an n-th 
non-linear Fourier-Pade approximant of (ax, . . . , a m ) with respect to ctq. 

Proof. In the proof we make use of multipoint Hermite-Pade approxi- 
mation. Fix n € Z™. For each j € {1, . . . , m}, choose an arbitrary set of 
| n | H — rij points contained in Ao 

Xn,j = {Xn,j,li ■ • • ; ^nj'^nl+nj ) ^ ^nj ; 
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where 

A n j = {(a?i, • • -,x\ nl+n .) G A n|+rtj : xi < • • • < x\ n \ +n .} . 

Let 

w n,j{x) = (x — £ n ,j,l) • ' ' (x — 37nJ,|n|4-n 3 ) j 

and consider the simultaneous multipoint Pade approximant which inter- 
polates the functions dj,j = l,...,m, at the zeros of w n j respective- 
ly. That is, On,i/<? /q n ) is a vector rational function such that 
deg(p n ,j) < |n| - 1, j = 1, ... ,m, deg(q n ) < |n|, q n 0, and 



(23) qri ° 3 Pn ' j G W(C \ supp((7,-)). 

From Lemma|2]we have (JJJ) and (jSJ) • Once we have determined q n , for each 
j = 1, . . . , m, we define the monic polynomial f2 nj , deg(finj) = l n l + n j, by 
the orthogonality relations 

(24) fy k n nA y) ( * - / <^^M^M ^| d(To(y) = o, 



ij(y)Qn,j(y) J y-x w nJ (x) 

k = 0, . . . , |n| + rij — 1. For each j = 1, . . . , m, these relations determine 
a unique O n j since the (varying) measures involved have constant sign on 
A . 

The polynomial 0, n j has exactly |n| + rij simple zeros in the interior of 
A . Set 

^nj = (yn,j,l7 • • • ; yn,j,\n\+rij ) £ ^n,j ; 

where y n j,i < • • • < y n ,j,|n|+rij are the zeros of ^ n ,j- 

Since for each j = 1, . . . ,m, the distance between Aj and Ao is greater 
than zero, the correspondence 

(-^n,l j • • • > -^n.m) * • • • j ^n,m) ) 

defines a continuous function from A nj i x • • • x A njjn into itself with the 
Euclidean norm. The continuity of this function is an easy consequence of 
the fact that Ao fl A,- = 0, j = 1, . . . ,m. By Brouwer's fixed point Theorem 
(see page 364 of 4 ) this function has at least one fixed point. Choose a 
fixed point. Then, w n j = ^ n j, j = l,...,m. Consequently (|2*i)) can be 
rewritten as 



(25) i^/^feiW*^)' 
■/ \<ii,j(y)<in,j(y) J y~ x w n,j{x) 



dcr (y) = 0. 



= 0, . . . , |n| + rij — 1, and taking into consideration (jHJ) we obtain that for 
each j = 1, . . . , m, 

f (a^x) - x k da j {x) = , k = 0, . . . , |n| + rij - 1 . 
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From the definition, it follows that (p n ,i/<7n, • • • >Pn m /<Zn) is an n th non linear 
Fourier-Pade approximant for the Angelesco system, taking S n j = Pn,j, 
j = 1, . . . , ra, and T n = q n . □ 

Let (J^j^i ■ ■ ■ i S ^p m be any non-linear Fourier-Pade approximant with 
respect to the Angelesco system (a±, . . . , a m ). From ii') it follows that 
<5j(z) — S t 3 (z) has at least |n| + rij sign changes on Ao- Let W n j be the 
monic polynomial whose zeros are the points where this function changes 
sign on Ao- Obviously, deg W n j > |n| + nj and 

(26) — e W(C\supp(cr i )), j = l,...,m, 

is analytic on the indicated region. (These polynomials W n j do not coincide 
with those of the linear case.) Using Lemma [2] it follows that 

(27) j x k ^^Lda 3 (x) = 0, fc = 0, l,..., nj -l, j = l,..., m. 
and 

where T nj - is the monic polynomial whose zeros are the rij zeros of T n lying 

in the interior of Aj. Combining (|28|) with ii') we obtain 

(29) 

J V T nj(y)|r n ,i(y)M \v-*\ \ W »M\ J 

The proof of Theorem |21 is similar to that of Theorem ^ First, we study 
the asymptotic zero distribution of the polynomials T n j and W n j. Then, 
we use this result to obtain the asymptotic behavior of the remainder in the 
approximation . 

Theorem 4. Let (do; 0"i, . . . , a m ) € Reg and consider the sequence of 
multi-indices A = A(pi, . . . ,p m ). Then, there exists a vector measure ~fi = 
(Tii, . . . ,~p2 m ) G -Mi such that for each j = 1, . . . ,m 

* H PJ v T n>j =~Pj, * lim vw . = ]I m+j . 

neA J neA J 

Moreover, Jl = ~p(C2) is the vector equilibrium measure determined by the 
matrix C2 on the system of intervals Fj = Aj, j = 1, . . . , m , Fj = Ao, j = 
m + 1, . . . , 2m . 

Proof. Let us show that the sequences of measures {z^T n ■ } and {v Wn }, n G 
A, have only one accumulation point. Let A' C A be a subsequence of indices 
such that for each j = 1, . . . , m 

* ne™ UTn ' j =Uj ' * ne™ UWn ' 3 = Um+j ' 
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(Notice that Vj G = 1, ...,m, and Uj G A4i(Ao),j = m + 

1, . . . , 2m.) Therefore, 

(30) lim \T ntj {z)\^ = exp(-y"i(z)) , 

neA' 

uniformly on compact subsets of C \ A j , and 



(31) lim |W n ,,-(z)| |n| +^ = exp(-V Um +Uz)), 

neA' 

uniformly on compact subsets of C \ Ao- 

As we have seen, T nj is orthogonal with respect to the varying measure 

l^jdcrj. Using and (JST}, we obtain 

lim — log ffi^M = .l+Wy^ (x) + y 2L V ^ (x) , 

uniformly in Aj . By (|27|) and Lemma^J is the unique equilibrium measure 
for the extremal problem 

(32) V"i (x) + ^ |^ fc (x) - ±±PlV Vm +i (x) > ^ , x G Aj , 
with equality for all x G supp(i/A and 



i 

2" , 



(33 » ^{l^^w 1 } 

These relations are completely similar to those obtained for the linear case 
(see p3p and (|14|)). On the other hand, W n j satisfies the orthogonality 
relations (|29[). We can apply once more Lemma |I] obtaining that, for each 
j = 1, . . . , m, v m+ j is the unique equilibrium measure for the extremal prob- 
lem 



(34) V Vm+i (x) - -?i—V v i (x) - -T^r^V* (x) > Vm+j , x G A 



with equality for all x G supp(f m +j). These relations differ from those ob- 
tained for the linear case (see (|17|)) 

If we look at the matrix corresponding to this system of equations we see 
that it is not symmetric. Let us rewrite the system as follows. Multiply 
equations (|3*2*)l times 2p 2 - and we obtain for each j = 1, . . . , m, 
(35) 

2pjV* (x)+^ PjPk V Uk (x)- Pj {l+ Pj )V Um +J(x) > 2-rijp) = Wj , x G Aj , 

k+j 

with equality for all x G supp(z^). With equations (|34"|) we have to work 
harder. First, let us multiply them times 2(1 + Pj) thus obtaining for each 
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j = 1, ... ,m, 

(36) -2 Pj V"i(x) - ^2p k V Vk {x) + 2(1 + Pj )V Vm +i(x) > 

k^j 

2r/ m+j (l +pj)= r]' m +j > x G A , 

with equality for all x € supp(^ m+J ). 

Let us show that in this second group of equations we have equality for 
all x G Ao- In fact, notice that 2(1 + pj)v m+ j is a measure on Ao of total 
mass equal to 2(1 + Pj)- On the other hand 

Mi 

is a measure of total mass pj + Yl^iPk = 1 + Pi < 2(1 + Pj) supported on 
the set U^L 1 Afc which is disjoint from Ao- Therefore, 

2(1 +pj)u m+j = (2pjVj + ^2pkVk)' + (1 +Pj)ua , 

where (•)' denotes the balayage onto Ao of the indicated measure and wa 
is the equilibrium measure on Ao (without external field), since these two 
measures are supported on all Ao so is their sum. Thus, supp(z/ m+ j) = Ao . 

The idea now is to take row transformations on the system of equations 
(|36|) to transform it conveniently. The matrix of this system of equations is 

/ -2pi -P2 ••• -p m 2(1 +pi) ••• \ 

- Pl -2 P2 ■■■ - Pm 2(1+ pa) ••• 

V -Pi -Pi ■■■ -2Pm ••• 2(1 +p m ) ) 

Since each column has a common factor we will carry out the operations 
without the common factor and afterwards place them back. Thus in 
columns k = 1, . . . ,m we factor out —pk and in columns k = m + 1, . . . , 2m 
we factor out 2(1 +pfc_ m ), respectively. The resulting matrix is 



/ 2 1 ••• 1 1 ••• 0\ 
1 2 ••• 1 1 ••• 

y 1 1 ••• 2 ••• 1 J 



B T ) 



where I denotes the identity matrix of order m. We know that the submatrix 
B is positive definite and through row operations can be reduced to the 
identity. This is the same as multiplying ( B X ) on the left by B^ 1 . 
Doing this we obtain the block matrix 



I B 



-i 
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It is easy to check that 



1 



m + 1 



m 
-1 



m 



-1 \ 
-1 



m 



J 



V-i -l • 

Multiplying back the factors we extracted we obtain the matrix 















-Pr, 



2m(l+pi) 


-2(l+p 2 ) 


-2(l+p m ) 


m+1 


m+1 


m+1 


-2(l+pi) 


2m(l+p 2 ) 


-2(l+p m ) 


m+1 


m+1 


m+1 


-2(1+Pl) 


-2(l+p 2 ) 


2m(l+p m ) 


m+1 


m+1 


m+1 



Therefore, the system of equations ()36|) is equivalent to 
(37) 



- Pj V^(x) + 2 "l H fV ^(x)- 



m + 1 



Mi 



where 

(^m+l ) • • • > ^m) = & (Vrn+l i 

Finally, multiply the jth equation in (|37|) times (1 + Pj) to obtain 



(38) 



-py(l + Pj)V^(x) + 2m i 1+ f j) V ^(x)- 



E 



2(l+p fc )(l+ Pj ) 
m + 1 



m+1 



V 1/m+ -» (x) = r/ m+j (1 + Pj ) = w m+j , x G A . 



The system of equilibrium problems defined by Q35|) and l|38|) has the inter- 
action matrix 

Cl,l C\2 



Co 



C 2j i c| 2 



defined in Section 1. Thus, the corresponding equilibrium problem has at 
least one solution given by (v%, . . . , u m ). According to Lemma 1, (v%,..., u m ) 
is uniquely determined if we prove that C2 is positive definite. 

Let us show that C2 is positive definite. The first m principal minors of 

(i) 

Ci and C2 coincide and we already know that they are positive. Let C 2 
denote the principal minor of C2 of order j where j £ {m + 1, . . . , 2m}. For 

(i) 

each k = 1, . . . , m, factor out from the kth row and fcth column of C\ . 
From the row and column m + k, k = 1, . . . , j — m, factor out 1 + In 
the resulting determinant, for each k = 1, . . . ,j — m, add the fcth row to the 
(m + /c)th row and then to the resulting determinant add the kth. column to 
the (m + k)th column. We obtain 

det(C { 2 j) ) = [ Pl ---p m (l+p 1 )---{l+p j ^ m )] 2 x 
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1 1 
1 1 
1 1 



1 1 



1 



1 


1 


1 


1 


1 


1 


1 

2m 


1 

m—1 


1 

m—1 


m+1 
m—1 


m+1 
2m 


m+1 
m—1 


m+1 


m+1 


m+1 


m—1 


m—1 


2m 


m+1 


m+1 


m+1 



In the determinant above, delete the row m + 1 from the following ones and 
in the resulting determinant add to the column m + 1 those after it and we 
get 



2 


1 


•• 1 




3 ~ 


- m 


1 


1 


1 


2 


•• 1 




j ~ 


- m 


1 


1 


1 


1 


•• 2 




j ~ 


- m 


1 


1 


1 



1 




•• 1 
•• 


(m+l)+(i 


— m)(m— 1) 


m—1 


m—1 




m+1 




m+1 
1 


m+1 










• • 










1 


2 1 




1 j — m 




2 1 •• 


■ 1 





1 2 




1 j — m 




1 2 •• 


• 1 













+ 








1 1 
1 1 




2 j — m 
1 j — m 




1 1 •• 
1 1 •• 


• 2 

• 1 




(m+l)-2(j- 

m+l 



(j — m) + (m + 1) — 2(j — m) = 2m + 1 — j > . 

With this we conclude the proof. □ 
We are ready to prove Theorem El 

Proof of Theorem 12 From (|28j) the asymptotic behavior of dj{z) — 

S ■ (z) 

T( z ) can be expressed in terms of that of the sequences of polynomials 
W^nj, T n j, and ( n j, where 



1 



mm 



\ T n,j{x) 

On account of Theorem we have 



Mx)? \^p\d*M , Q{x) 



2 \ T n,j(x)\d<Jj(x) 



\WnAx)\ 



(39) 



lim |Wnj(^) 



VIH 



neA 



exp{-(l +p j )V^™+j(z)} , 
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uniformly on compact subsets of C \ Ao, and 

(40) hm l^-OOp/H = exp{-2 Pj V^(z)} , 

uniformly on compact subsets of C \ Aj , where Jl = 7* (£2)- Using ()33[). we 
have 

/ 1 y/M 

(41) hm \ 7T~ J = ex P{- 2 Pj r lj} = ^P{-Wj/pj}. 

Combining (j2SJ), ([HU]). (0DJ), and flUJ), we obtain that (J3J) holds true uniformly 
on compact subsets of the indicated region. □ 

4. Comments on Lemma [T] 
Let Ai(Fk), k = 1, . . . ,N, be the class of all measures on Ai(Ff.) and 
M = M ( F 1 ) x ••• x^(Fjv). 
Define the mutual energy of two vector measures (J- 1 ,^ 2 by 

N 

(42) J( M 1 ,M 2 )=E / I c^lny^—d^zWkix). 

j,k=l J J \ z xl 

The energy of the vector measure /i € is 

v 

(43) J(aO = ^2 c j,kI{^j,Vk) , 

i,fc=i 

where 



(44) I(fij,fi k ) = J J In * dixj(z)dn k {x) ■ 

Given a real symmetric positive definite matrix C and fj, E M define the 
combined potentials W^j = 1, . . . , N, as in the introduction, and the vector 
potential = (W^ , . . . , WJy). These formulas may be rewritten as 



(45) J(mV) = J W* (zW{z), 

where 

| W»\zW l (z)=jrj W?\z)dn}(z), 



and 

(46) J(fi) = J W fl (z)dn(z) . 

If fj,, fj, 1 , (J? £ £ are vector charges whose components have finite energy, the 
energy of a charge J(/u) and the mutual energies of the charges J(/i 1 ,/x 2 ) 
can be defined analogously by formulas (|4^|) and (|42|) . respectively. 

If c j>k > when Fj D F k ^ $,j,k G {1, . . . , N} the functional J^ 1 , /x 2 ) 
and J(/i) are lower semicontinuous in the weak topology M. (see Proposition 



20 M. BELLO-HERNANDEZ, G. LOPEZ-LAGIMASINO, AND J. MlNGUEZ-CENICEROS 



5.4.1 in [3]). Consequently, the functional </(//) attains its minimum in Mi. 
In Proposition 5.4.2, using a unitary decomposition of C, the authors prove 
that J(/u) is a nonsingular positive definite quadratic form on the linear 
space 8. In Proposition 5.4.2, the extra condition on the coefficients of C 
is not needed. Therefore, we can say that if there is a minimizing vector 
measure then it is unique. 

Let < e < 1 and ii , /i 2 S Mi. Assume that the components of fi 1 ,^ 2 
have finite energy. Set fx = efi 2 + (1 — e)^ 1 G Mi. It is algebraically 
straightforward to verify that 

(47) J(jZ) - J(/i 1 ) = e 2 J(/x 2 - /i 1 ) + 2e f {x)d(iL 2 - ^){x) . 



Assume that J(/u , ) minimizes the energy functional. Dividing by e and 
letting e tend to zero, it follows that 

(48) f W^{x)d{n 2 - n l )(x) > 0. 



for all /U 2 € M.\. Reciprocally, assume that (|l8]) takes place for all \i 2 € Mi, 
then using (|47|1 it follows that [i 1 minimizes the energy functional since 
J(^ 2 ~ M 1 ) > for all fj, 1 ,^ 2 G C. 

Now, let /Z € A4i be a solution of the equilibrium potential problem 
determined by C on the system of intervals Fj ,j = 1, . . . ,N . That is 

Wf(x) =wf, x G supp(7i j ) , 

where Wj = mi{Wj{x) : x £ Fj}. Hence, for all /i € 

AT JV 

W"(x)ci0x-7Z)( ! r) = X; / ^(xK^--7Z ; -)(»)>Z) u '?- u f 
i=i 7 i=i 



and it follows that /i minimizes the energy functional. With this we conclude 
the comments on Lemma 1. 
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